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By Max. A. Heaslet and Franklyn B. Fuller 


SUMMARY 

The minimization of inviscul fluid drag is studied for aero- 
dynamic shapes satisfying the conditions of linearized theory , 
and subject to imposed constraints on lift, pitching moment , 
base area , or volume. The problem is transformed to one of 
determining two-dimensional potential flows satisfying either 
Laplace's or Poisson's equations with boundary values fixed by 
the imposed conditions. A general method for determining 
integral relations between perturbation velocity components is 
developed. This analysis is not restricted in application to 
optimum cases; it may be used for any supersonic wing problem. 

For given base area , general formulas are found that cover as 
special cases quasi-cylindrical bodies of revolution , wings 
having plan forms with fore-and-aft symmetry , slender bodies , 
and certain classes of yawed wings. The drag can in fact be 
determined from a unidimensional flow analysis in a duct of 
known shape. For given volume , minimum wave drag of a 
ducted body of revolution of arbitrary radius is expressed in 
closed analytic form. The elliptic wing is treated, and a pos- 
sible source of difficulty connected with unreal shapes in given 
volume problems is found. In the case of ducted bodies of revo- 
lution, the singularity distribution corresponding to the mini- 
mum drag is determined. 

Particularly simple results are found for a family of wings 
with curved leading edges with lift specified and center of pres- 
sure fixed at the 60-percent-cliord position. General relations 
for integrated loading along oblique lines are derived for this 
family of wings. 

INTRODUCTION 

To seek conditions under which the wave and vortex drag 
of a given wing or body is minimized is to seek conditions 
for economical supersonic flight. It is also a common ex- 
perience, in the study of such problems, to find that a gratu- 
itous economy appears to affect the analysis itself. Almost 
invariably, simplicity characterizes the final forms of the 
results in comparison with predictions carried out for wings 
and bodies chosen with less discrimination. In the present 
paper, the minimization of wave and vortex drag for various 
aerodynamic shapes is studied. Some side conditions, such 
as given lift, given volume, etc., must be specified in order 
to set a meaningful variational problem for the shape at 
hand. The side condition of given base area is noteworthy 
as leading to results with the simplicity mentioned above, 
for, as will be shown, the general expression for minimum 
drag assumes the most elementary form possible while at 
the same time retaining the relevant parameters and being 
dimensionally correct. 


The starting point of the present work was the expression 
for drag given by G. N. Ward (ref. 1) in his study of thin 
lifting bodies, that is, wings and bodies for which linearized 
supersonic flow theory applies. Work of this type was first 
reported by Nicolsk}^ in reference 2, and detailed results 
(refs. 3 and 4) have since become available. The body 
shape is assumed to be enclosed b}- a characteristic surface 
generated as the envelope of both the downstream-facing 
Mach cones, with vertices on the forward edge of the body, 
and the upstream-facing Mach cones, with vertices on the 
trailing edge of the body. The drag (wave plus vortex) 
is then given by a surface integral of the induced velocities 
over the downstream portion of the Mach envelope, and 
other forces, moments, and geometrical properties are 
similarly determined. This particular control surface has 
analytical advantages similar to those exploited by II. T. 
Jones (refs. 5 and 6) in the use of combined flow" fields. 
Jones adopts a perturbation potential equal to the sum of 
the potentials in forward and reverse flow. He then shows, 
for example, that for a plan form of given base area the 
necessary condition for minimum wave drag is that the 
pressure in the combined flow field be a constant over the 
plan form. It follows that locally the combined flow 
potential is a two-dimensional harmonic function. The 
direct use of combined flow fields has been further extended 
and applied by Graham (ref. 7), and is used in combination 
with the control surface mentioned above by Germain 
(refs. 8 and 9). Along the Mach envelope used by Ward 
the perturbation potential in forward flow- is equal to its 
value in the combined flow field. Drag minimization then 
determines conditions on the control surface, and the 
potential on this surface differs from a two-dimensional 
harmonic function by a known amount. The conventional 
perturbation potential is retained, but the determination 
of the body shape is still not direct. Mathematically, one 
needs to invert an integral equation, and questions as to 
existence and uniqueness of the solution arise. In the case 
of ducted bodies of revolution of given base area or volume, 
the integral equation will be solved explicity in a later 
section. 

When the expressions for the forces and moments and 
geometrical quantities in terms of integrals over the same 
control surface are known, it is possible to combine the drag 
with one or more of the others, as constraints, and set up 
an optimization problem. Variational methods then yield 
the result that, in general, the potential in the control surface 
must satisfy a Poisson equation in the lateral coordinates. 
Inspection leads to a number of cases in which solution by 
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analytical methods is possible, and recourse to numerical 
techniques is not required, though it is of course available 
for more difficult problems. 

Finally, a method for generating integral relations among 
perturbation velocity components is outlined. By these, it 
is possible to gain knowledge of aerodynamic loading on the 
wing integrated along oblique lines. So long as the obliquity 
of the lines is such that the component of stream velocity 
normal to them is supersonic, the variation of the integrated 
loading in the stream direction is obtainable. The so-called 
chord loading is a special case of these results. The span 
loading is a by-product of the solution of the variational 
problem, for this involves determination of the perturbation 
potential in the rear Mach surface springing from the 
trailing edge. Then, when the perturbation potential on 
the trailing edge is known, the span loading is found directly. 
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IMPORTANT SYMBOLS 

base area of wing or body 
drag coefficient 
lift coefficient 

curves bounding a region in an x= const, plane (see 
fig. 2) 

drag of wing or body in a supersonic flow 
function defining (by x=J{y 1 z)) a characteristic sur- 
face springing from a trailing edge 
streamwise extent of wing or body 
lift 

integral of local loading f — ■ j along oblique line (see 
eq. (66)) 

tangent of angle of sweep or yaw 
Mach number in the free stream 
pitching moment, positive for a nose-up moment, 
taken about the line x=x m , y=z = 0 
inner normal to a plane curve 
inner normal to a surface 
pressure 

load coefficient (upper-surface pressure minus 
lower-surface pressure divided by free-stream 
dynamic pressure) 
local velocity 

free-stream dynamic pressure, ^ p m UJ 

radial coordinate in axially-symmetric problems 
arc length 

region in an x= const, plane (see fig. 2) 
perturbation velocity components in x,y,z direc- 
tions 

speed of free stream along x axis 
volume of a wing or body 
Cartesian coordinates 
MJ-l 

space curve (defined in fig. 1) 

Lagrange multipliers (see eq. (11)) 
conormal (see eq. (18a)) 
density in the free stream 

^surfaces 


(p perturbation velocity potential 

X perturbation potential on surface x=f(y,z ); x= 

<p[f(y,z),y,z\ 

d 2 d 2 

V“ two-dimensional Laplace operator, 

1 ox 2 oy 2 

ANALYSIS 

The following analysis is divided into several sections. 
First, the geometrical quantities, base area and volume, are 
expressed in terms of integrals over a given control surface. 
Then relations for the lift, drag, and pitching moment are 
given. These results can be combined in various ways to 
set up variational problems where drag is minimized while 
one or more of the remaining quantities is held fixed. 
Finally, a method by which integral relations among the 
perturbation components can be derived is outlined. 

In all this work, it is assumed that supersonic small- 
disturbance theory applies. There is then a perturbation 
potential function <p(x,y,z) satisfying the equation 

P 2 <Prr—<Pvu— i Pzz=Q ( 1 ) 

where 


p*=MJ - 1 

The velocity components in the coordinate directions are 
U m -\- (p x ,ipy,Vz- It will be assumed that the given body can 
be represented by a planar or cylindrical reference surface, 
the latter having directrices parallel to Ox, the free-stream 
direction. Boundary conditions are to be satisfied on the 
reference surface. In order to avoid difficulties concerning 
gaps or holes in the body surface, it will be assumed that 
unique leading and trailing edges exist, and that the thick- 
ness distribution function does not vanish between these 
extremities. 

RELATIONS FOR BASE AREA AND VOLUME 


Consider a cylindrical reference surface as mentioned 
above, and draw in the enclosing characteristic surfaces as 
shown in figure 1. The reference surface itself is denoted 
2 0 . The envelope of Mach cones springing from the lead- 
ing edge is 2,, and the envelope of Mach cones from the 
trailing edge is S 2 . These surfaces, and 2 2 , intersect 
along a space curve r 2 , as shown in figure 1 . The relations 


z 



Figure 1. — "Reference surface and characteristic enveloping surfaces. 
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for base area and volume of the body represented by the 
reference surface 2 0 will be derived by applying Green's 
theorem to the volume bounded by the surfaces 2 0 , 2i, 
and 2 2 . 

The following form of Green's theorem will be used: 


(It 


fff F(yU ) • (VF) dr=-f f FV^dZ-fff F[V 

T 

2 o-}- 2 i-f-Z 2 T 

whore 


T 

volume bounded by 2 0 +2! + 2 2 

F, U, V 

arbitrary scalar functions of position 

N 

interior normal to a surface 

V 

Now sot 

vector differential operator, grad 


F= P =p m (\—MJ to first order 

U=$=U m x+<p(x } y,z) = total potential 
V=x k 

The integral relation becomes 

JJJ 'l>(U.+p x )kt*- 1 dr= -JJ, px ‘ * ^ds -fff >[V- (pV*)] 
V 2 o+Si + 22 t 


dr 


and the last integral on the right vanishes by virtue of the 
continuity equation for steady flow. The remaining quan- 
tities are 

JJJ***-‘p(tf„+«,) <*r=- JJpx* ~ dZ (2a) 


So+S, + S 2 


Take first the case &=(). The left side of equation (2a) 
vanishes and one has 


if'iwjj'Mf 

So 2i S 2 

The integral over 2 0 is, to first order, 
d4> 


o < 2b ) 


If* 


d;V 


dZ=pJU 


Jfr ^ p r JJ ' /v 


where Z=Z(x,y ) is the equation of the surface of the wing or 
body. The integral therefore gives the increment in frontal 
area of the wing or body between the leading and trailing 
edges, symbolized here as A; 


A 


=JJ 

So 


* dZ 
dx 


dx dy 


( 3 ) 


Next, consider the integrals over 2 X and 2 2 . Let the surface 
2 2 be defined by the equation 

x=f(y, 2 ) 

The normal N to 2 2 has direction cosines 

M • M • AT 1 . fv . f* 


and the normal to 2 X has clearly Ni=l/M n . Thus equation 
(2b) becomes (since <p= 0 on 20 

p JJ„A=p a ff (.—P 2 <Pz—fvVu—fz<Pz)dy dz 


An essential simplification of this formula results if one 
introduces the function x(y>z) where 

x=<p\f(y,z),y> 2 \ 

and x=f(y,z) is the equation of the rear characteristic 
surface 2 2 . Then, since 

Xv= <Pxfy+ <Pv J Xz= <Pxfz+ <Pz 

and jy + fz 2 =P 2 , the relation for base area becomes 

A=—~ jj {fyXy+f^dydz (4a) 

& 

In equation (4a), S is the projection, shown in figure 2, of 
the surface 2 2 on an x= const, plane. As shown also, Ci is 



plane. 

its outer boundary, or the projection of the space curve T t 
shown in figure 1, and C 2 is the projection of the reference 
surface of the wing or body. By application of Green's 
theorem, and the knowledge that x=0 on T u hence on C\, 
equation (4a) can be rewritten in the forms 


i* df . 

- |x 5i* 


A ~u~Jj xv ' fdydz+ irJ 

s c 2 

-rhi! fv ' xdyd ’ + rj f ^ d ‘ 


(4b) 

(4c) 


Ci+ft 
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where 

n interior normal to curve in x— const, plane 
s arc length in x= const, plane 

The formula relating the volume of a wing or body to 
integrals of the function x follows by taking k 1 in equation 
(2a) ; one has 

V =J fh(f,s)ds-~jj x(v 2 P+ 2 p)dy dz-^j- j X <h 

Cl »S' C2 

( 5 ) 

where 

1" incremental volume of wing or body 

h(fjS) deviation of body from control surface at trailing 
edge 

Equation (5) gives the volume of bodies that are either open 
or closed at the base. The first integral on the right of 
equation (5) may vanish for several reasons aside from the 
obvious case of zero base area; in these cases the volume is 
expressed as 

v —m If f a * -£ * «*> 

5 


s s 



It is worth noting that since neither base area nor volume 
can be influenced by position of the center of coordinates, it 
is possible sometimes to effect certain economies in algebraic 
manipulation by a judicious choice of origin. 

FORCE AND MOMENT RELATIONS 

Expressions for lift and drag forces are obtained by appli- 
cation of the law of conservation of momentum in the 
volume bounded by 2 0 , and Z 2 (see fig. 1). The vector 
force relation is 

-> rr -> c 

pNdZ— JJ pq(q-N)dX 

2i+2o Si+Sj 

where 


J'=—pJ\jjx s </ydz (7a) 

.S’ 

= l T » J X cos (z,ri) ds (7b) 

z> =— yJJ XV t xdydz-^J x * x -ds (8a) 

s 2 

= " 2 "J] (tf+x*)dydz (8b) 

s 

The pitching moment is next found. The vector moment 
relation is 

M — £| p(q-N)(r m Xq)d2— JJ p(r m XA0d2 

2,+S. 2,+Z 2 

where r m is the vector distance between the moment center 
and an integration element on the control surface. If only 
the pitching moment is considered, with nose-up pitching 
moment taken as positive, and moment center at (x m , 0, 0), 
the linearized result is 

M = — pJJ* J | XV 2 (zJ)dydz— 



If lift L= 0 or r m =0 so that pitching moment is calculated 
about the origin, one finds 

M =-p„U„ffxv*(zf)dydz-pM„j x d( £ ] ds (10a) 

s C2 

= — zjv 2 xdydz— pJJ„\j zj ^ ds (10b) 

VARIATIONAL PROBLEMS 

The problem of minimizing drag under the constraint of 
given base area, or volume, or lift, or pitching moment can be 
set up with the aid of equations (4), (6), (7), (8), and (10). 
It is only necessary to apply standard variational procedure 
to any of the expressions 


F vector force on body 
-> 

q local velocity vector 

— ) 

N unit inner normal to surface 

p local pressure 

p local density 

The force can be resolved into lift and drag components, and 
the results consistent with linearized theory are (as given in 
ref. 1) 


I^D-XA" 

i 2 =d+»v ^ 

I z =D—aL ^ 
I±—D — tM ^ 


( 11 ) 


where X, p, a, t are Lagrange multipliers. In the language of 
combined flow fields (refs. 5 and 6) the Lagrange multipliers 
yielding minimum drag can be identified with constant values 
of the longitudinal and vertical velocity components or their 
gradients. In this way, for example, X and p are, respectively, 
constant pressure and pressure gradient in the combined 
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field. Less general interpretations may possibly be given the 
Lagrange multipliers in specific classes of problems. One 
such interpretation is given in the discussion following 
equation (40) of this report. 

It is possible to combine the variational problems with 
each other. For example, if it be required to find minimum 
drag with given base area and lift, the quantity to be mini- 
mized would be D — \A — aL and so forth. The results found 
by applying the variational procedure to equations (11) will 
be given next. Each is a two-dimensional flow problem in the 
lateral variables y,z. 

Given base area: 


d 

dn 


+ 

■O 1 

8 cL 
8 1 

II 

O 

in S 

8 

II 

O 

on C 2 

x=0 

on Ci 


d=\a 


Given volume ( zero base area): 

v ’[ x+ i 

)J=0 on C 2 


i=° 


X=0 


on C Xj 


(12a) 


(12b) 


(13a) 




(13b) 


upon the boundary curve C\. Several cases where exact 
solutions are obtainable will be treated in later sections. 
However, recourse to approximate methods is usually indi- 
cated. One approximation in which the wing is distorted 
slightly in order to obtain a boundary curve C\ for which the 
two-dimensional problem is solvable is discussed in refer- 
ence 10. Germain (ref. 8) has used development in series. 

INTEGRAL RELATIONS BETWEEN PERTURBATION 
VELOCITY COMPONENTS 

Some useful and interesting relations between integrals of 
the potential function taken across the reference surface 
2 0 to integrals of x in S 2 (see fig. 1) will now be derived. 
These results are not specialized to minimum-drag solutions, 
but apply to any case. Consider a planar, supersonic- 
edged wing, and choose a line in the wing plane (z= 0) that 
cuts the x axis in (z 0 ,0,0) and which makes an angle y 0 with 
the y axis; 

x—y tan ii 0 =x 0 (16) 

Further, let the angle no be such that tan making the 

normal component of free-stream velocity always supersonic. 
This line will always lie in the supersonic “zone of silence” 
corresponding to the point (x 0) 0, 0). Under this condition, 
a plane, S 3 , tangent to the Mach cone springing from 
(x 0 , 0, 0) can be made to pass through the line given by 
equation (16). This plane is (in a notation introduced by 
Hayes, ref. 11) 

x— ft cos dy-\-& sin 6z=x 0 (17) 

where 13 cos 0=tan n 0 . The resulting situation is shown in 
figure 3 for the case of a sonic-edged wing. 


where k is a constant to be determined by application of the 
zero base area condition. 


Given lift: 


V 2 X=0 in aS° 
(x+f7„crz)= 0 Oil C 2 - 
X=0 on 


(14a) 



(14b) 


Given pitching moment ( zero lift) : 

V 2 (x+C/ a> rz/) = 0 in S' 

x=0 o\\C x> 


(15a) 


-y 



Figure 3. — Characteristic plane cutting wing and enveloping surfaces. 


D=-^ M (15b) 

where k is a constant to be determined by the condition of 
zero lift. 

In each of the problems listed in equations (12) through 
(15), the possibility of obtaining a solution depends largely 


If Greeks theorem for equation (1) is applied to an 
arbitrary region enclosed by a continuous surface 2, the 
fundamental integral relation 


(18a) 
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results. The derivative dip/dv is the gradient of potential 
along the conormal v with direction cosines v u v 2 , v 3 that are 
related to the direction cosines n h n 2 , n 3 of the inner normal 
to the surface 2 by the equations 

—P 2 ni= n 2 = \v 2 , n 3 = \v 3 

Consider now the region bounded by 2 0 , 2 U and S 3 , and 
in which z ^> 0 (in fig. 3, 2 0 , 2 2 , are defined as in fig. 1; 

2 3 is the plane of equation (17)). On 2i and 2 3 , it is easily 
found that A = 0, and on 2 0 , A=1 For this region, equation 
(18a) becomes 


SJJM/tfHT 

2o 2.) 2,3 


d<p 

dv 


( 12 = 0 


2o 23 


^r/2=0 

dv 


i jj !>, 

2: 0 


(18b) 


Since the conormal v lies in the surface 2 U (and also in 2 3 ), 
the integral over Si will vanish because ip is zero thereon. 
Further, on 2 0 , d/dv=d/dz and d2 = dx dy. Equation (18b) 
now becomes 


(19a) 


The integral over S 3 in equation (19a) will be greatly 
simplified if one takes d2=dv r/.s* , where ds is an element of 
length normal to dv , and lying in the plane S 3 . It can be 
shown that under these conditions, ds always lies in an 
x= const, plane. Equation (19a) now' becomes 


(19b) 


where the line CD is the intersection of the plane given by 
equation (17) with the wing plan form. 

Next, one can repeat this analysis for the region ahead 
of S 3 , within 2, and S 2 , and below the wing plane. This 
time, the surface 2 2 will form a portion of the bounding 
surface, and a line integral across S 2 will result. One can then 
combine the relations obtained for the two regions and 
relate an integral of potential across the wing surface to one 
across the surface 2 2 . A detailed application will be made 
in a later section. 

APPLICATIONS 

The analysis discussed in the previous sections will now 
be applied to some particular problems. These will include 
a number of cases in which the configuration of the wing or 
body permits simple analytic solution of one or more of the 
variational problems listed as equations (12) through (15). 
Quasi-cylindrical bodies of revolution were mentioned in the 
Introduction as constituting a class of shapes for which a 
full solution to the minimum-drag problem is available. 
The simulating axial source distribution, from which surface 
shape can be calculated, is derived herein, leading, inci- 
dentally, to an interesting identity involving integrals of 
elliptic integrals. Finally, the results pertaining to integrals 
of wing loading along oblique lines across the wing are 
applied to a particular family of minimum-drag wings. 


WINGS AND BODIES WITH GIVEN BASE AREA 

By combining equations (8a) and (12a) and using Green's 
formula, one gets for the drag of an optimum wing or body 
with the given base area 

"-£// W+/,*^=V- X/(s+^ 50* 

s 1 

0 2 X 2 x r ,/dx , X df\ , 

-4^ s+ 2r. ..),/( (20) 

where q is the free-stream dynamic pressure, (l/2)p 00 ?7 00 2 . 
Further evaluation of the drag by equation (20) requires 
explicit knowledge of the function x, obtained by solving the 
potential problem of equations (12a). A large and particu- 
larly interesting class of wings and bodies for which the 
solution is immediate is characterized by the condition that 
/=const. on C x . This implies that the outer rim T u the 
intersection of the characteristic envelopes in figure 1, lies 
in a plane normal to the free-stream direction. For example, 
all wings with plan forms having fore-and-aft symmetry 
satisfy this requirement as do also all pointed configurations 
with subsonic edges so long as the nose and tail vertices 
determine a line parallel to the free-stream direction. 

In such cases, the solution of equations (12a) is 

X = 0 

and equation (20) gives for the drag 


D=^ S 
4<Z» 


( 21 ) 


By combining this with equation (12b) and eliminating X, 
one gets 

I) A 2 

~(^~pS (22) 

The simplicity of equation (22) is remarkable, and ex- 
amples of its diverse applicability are given below. Before 
proceeding to these applications however, it should be noted 
that a similar result applies to all planar wings whose en- 
veloping surfaces and 2 2 intersect in any plane parallel 
to the 3 axis. In this case, 

x=f=m(y+b 0 ) 

on the curve (\, where m«0) is the slope of the plane of 
Pi relative to the stream direction. The solution of equa- 
tions (12a) once more gives and then, for the drag 

from equation (20), 


/mV C df (3 2 \ 2 S__\ 2 m 2 S 

4g» Jc, J bn US 4 h/ 




Again eliminating X with the aid of equation (12b), one has 
D___ A . 2 


(24) 
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Ducted body of revolution with prescribed end diame- 
ters. This problem lias been considered by Parker in refer- 
ence 12, and by Nicolsky (see ref. 2). As shown in figure 4, 


r 



Figure 4. — Optimum body with prescribed end radii. 

a shape with minimum external wave drag is constructed so 
as to have an initial radius R 2 and a final radius R x . In order 
that the previous linear theory should apply, the restriction 
is made that the ratio 1 7?j — /? 2 1 should be a small quantity. 
If the origin of axes is in the front face of the body, the 
fore-and-aft Mach surfaces are 


x~i Q(r—R 2 ), x—l= — P(r— Ri) 
and the curve Ci is a circle of radius R 0 where 

Ro= (Z+/3/?i+/3i? 2 )/2/3 

From equation (22) drag is 

D 4(/A 2 -i?2 2 )V 

(Zoo (l+pR 1 +pR 2 ) 2 -^ 2 Ri 2 


(25a) 


Equation (25a) is of particular interest since it represents a 
whole spectrum of results that extends from slender-body 
theory, for pRJl and pR 2 /l small, to two-dimensional theory, 
for pRJl and @R 2 /l large. The slender-body result leads 
directly to the familiar Karman ogive formula (ref. 13), 


Dk = 4A* 
(Z co 


(25b) 


The two-dimensional result applies to the upper half of an 
optimum wing (a flat plate), and is found to be 


^ D 2a 2 
( D torBlq. p 


(25c.) 


where R is mean radius and a=(R 1 —R 2 )/l. 


Elliptic plan form with afterbody. — The problem of given 
base area along the rear edge of an elliptic wing was consid- 
ered first by R. T. Jones (ref. 6). The figure is a semi-infinite 
body with a cylindrical shape drawn downstream of the 

483644—59 2 


rear edge of an ellipse, see figure 5. 
plan form outline is assumed to be 


cr 6 2 


The equation of the 


z 



and the enveloping Mach surfaces are determined completely 
by the fore-and-aft Macli cones with vertices along the 

supersonic-edged portion of the plan form ^ where 
that is, the abscissas of the vertices lie within the region 
\x\ <a 2 /(a 2 +b 2 p 2 )y> 

The curve C\ has the equation 

?/ 2 z 2 

jL _l 


[(a*+b 2 0 2 )m 2 W) 2 


and is an ellipse with foci at (±6,0). Equation (22) then 
yields 

D A 2 . , 

2. m{a 2 +b 2 p 2 )» j 

If drag coefficient C D is based on plan-form area, equation 
(26a) can be re-expressed as 

Y 




7 [P 2 + (4/ttAR) 


?^21 H 


(26b) 


since the aspect ratio of the wing is JR= (46) /(xa). Perhaps 
the most convenient formula for comparison follows from 
equations (25b) and (26a) if the drag of the wing is expressed 
in terms of the drag of a Karman ogive with the same length 
and base area. The ratio is given by 

Da 1 


D k (a 2 -\~b 2 p 2 ) yi [1 + (ttPJR/4:) 2 ]* 


(27) 


The wave drag of the elliptic wing with cylindrical after- 
body, in the limit as aspect ratio approaches zero, is equal 
to the drag of the Karman ogive and afterbody. For finite 
values of aspect ratio, the wave drag of the flat wing is 
smaller than that of the body of revolution, the initial devi- 
ation of the ratio from unity being proportional to (pdR) 2 . 

Tapered plan form with afterbody. — As a third example, 
consider a plan form of arbitrary taper ratio with base area 
along its trailing edge fixed. When the trailing edge is of 
subsonic type a cylindrical afterbody is assumed added. 
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As shown in figure 6, root chord is equal to 2 a and span 
equal to 26. The tip chord is 2 d so that taper ratio X Q and 
aspect ratio At may be introduced in the form 

\ 0 =d/a 

Al=2b/[a(l + \ 0 )} 

So long as the leading edge of the wing is supersonic the 
characteristic trace C\ is as shown in figure 7 and is composed 
of arcs of circles and straight lines, the radii of the inner and 
outer circles being a/p and d/p, respectively. The distance 
between the centers of the two outer circles is 26. Once the 
leading edge of the plan form is subsonic, the central circle 
of figure 7 blankets the other parts of the figure and C x is 
the circle of radius a/p. 



Figure 7.— Characteristic trace for wing with tapered plan form. 



The area S is the sum of elementary geometric areas and 
is given by 


>= 2a 2 Hr 

* P 2 L-2 


v+ a+A °) 2 cog a+a(1 _ V) j 


where a, shown in figure 7, is given by 

r 


a =* 


arc sin 2(1— X„) ^ 0M( 1+X„) 

9 > 2(1— X 0 ) §:/3JR(l+Xo) 


From equations (22) and 25b), the minimum drag of the 
tapered wing relative to the drag of the Karman ogive of 
equal length and base area is 


7T 


° K t +( 1+\,)[(1 +X 0 ) 2 /? 2 J? 2 -4 (1— Xo) 2 ]^— 2(1 -V) arc cos 

when 

2(l-Xo)^^(l+Xo) 

^=1 

D k 

when 

2(1-X„)£/3JV(1+X„) 


(28) 


Special cases of interest are : 
Rectangular wing (X 0 =l) 


D _ 1 

T) k 1 H- (4/3 JR) /t 


(29) 


Diamond wing (X 0 =0) 


D 

D k 


l+I[(/3^) 2 -4] 

7 T 


1 



7T 





2 ^ pJSR 


r 


(30) 


= 1 , PMI£ 2 

Figure 8 shows a plot of D/D K against PAR for the elliptic, 
rectangular, tapered, and diamond plan forms. Base area 


and length of the wings are equal to these parameters for 
the Karman ogive. For large values of PAR the relative 
drag decreases as 1 /PAR. As the wing plan forms become 
slender, drag of the elliptic and rectangular wings approaches 
in the limit of vanishing pAl the drag of the ogive. The 
tapered wing, on the other hand, has a value of drag equal 
to that of the ogive for all values of taper ratio and aspect 
ratio satisfying the inequality 2(1 — X 0 ) ^pAR (1 + X 0 ). This 
relation is satisfied so long as the edges of the wing are sub- 
sonic. Changes in sweep angle of the leading and trailing 
edges produce no further change in the minimum drag of 
the configuration so long as the base area is held fixed. The 
value D k is the minimum drag for all such configurations 
lying within the fore-and-aft Mach cones from the nose and 
tail of the wing. Except in the case of the rectangular 
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plan form, the curves of DID K have zero slope at their 
peak values. 

Yawed elliptic plan form with afterbody. — For given base 
area, the drag of an elliptic plan form at angle of yaw xp can 
be calculated from equation (24). In order to justify this 
statement it is sufficient to show that the characteristic 
curve T x lies in the plane x=f=m(y J rb 0 ). The trace of 
r 1 ! in a yz plane is, in fact, another ellipse and the dimensional 
relationships between the plan form and the trace are as 
shown in figure 9. It is convenient in the derivation of 
these results to proceed inversely and to determine the 
plan form as an envelope of curves given by the intersections 
in the xy plane of fore-and-aft facing cones with vertices 
on IV Since the stream wise position of the origin is of no 
direct significance, the space curve r x may be assumed 
given bv 

B * C \ (31) 

x=my, J 

The Mach cones with vertices at the point (x u y u z x ) on T x are 
(x—x 1 ) 2 =P 2 [(y—y 1 ) 2 + (z— z x ) 2 ] (32) 

where 

Xi=my u 

The parametric equation of the envelope in the yz plane 
is found from equation (32) with z=0 and the y x derivative 
of the same expression, that is, by elimination of y x from the 
relations 

B 2 [(x— my x ) 2 — p 2 (y— y l ) 2 ]=P 2 C 2 {B 2 — y x 2 ) 
[B 2 m 2 -p 2 (B 2 -0 2 )]y x =(mx-p 2 y)B 2 
The envelope is, therefore, 


Equation (33) represents an ellipse so Jong as the initially 
chosen m satisfies the inequality 

m 2 <p 2 (B 2 -C 2 )/B 2 (34) 

The elliptical plan form is fixed by its major and minor 
axes and angle of yaw. The relationship between the plan 
form and trace curve is more conveniently carried out, 
however, in terms of the three quantities Z, b, 8 where l is 
streamwise length of the plan form, 2 b its width, and x=8y 
is the line passing through the points on the plan form 
where y= ±6. 

Elementary calculations performed with equation (33) 
yield the following relations 


l/2=(B 2 m 2 +l3 2 C 2 )» 

(35a) 

b={B 2 - -C 2 )* 

(35b) 

mB 2 
5 B 2 —C 2 

(35c) 

— 2 mB 2 

(B 2 -C 2 )-{m 2 B 2 +0 2 C 2 ) 

(35d) 


Zb — 



h 2 B H 


(B 2 -C 2 )x 2 -2mB 2 xy+ (m 2 B 2 +0 2 C 2 )y 2 = C 2 {0 2 (B 2 - C 2 )-m 2 B 2 ] 

(33) 


Figure 9. — Yawed elliptic wing and characteristic trace. 
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Iii figure 9 the plan form is also circumscribed by a 
parallelogram with sides inclined at the Mach angle. The 
equations of these lines are 

x=i Sy± (p—m)B, x= ~Py± ( /3+m)B 

from which it follows that their outermost intersection points 
are at y=±B and the line connecting the intersection points 
is x—my. 

The above results thus show that the Mach lines circum- 
scribing the plan form can be used to determine the span of 
the trace of Tj and the angle of inclination of the plane of 
IV The span of the plan form is, moreover, equal to the 
distance between the foci of the elliptic trace. 

From equations (24) and (25b) the drag of the wing and 
afterbody relative to the drag of the Karman ogive of equal 
length and base area is given by 

D l 2 _B*m 2 +i3 2 C 2 

D k 4BC(P 2 -m 2 ) BC(JP—m 2 ) 

The results can be summarized as follows 

D f l+jS 2 f+[(l+^ 2 ) 2 — 4 0 2 5 2 £ 4 PO » 1 

W K ~~ \ 1 -0 2 £ 2 + [ (1 + 0 2 £ 2 ) 2 — 4/3 2 5 2 £ 4 1 ' J [ (1 +/3 2 £ 2 ) 2 — 4/3 2 5 2 £ 4 ] * 

(37) 


equation (37) are of interest in showing the variation of drag. 
For example, when 5=0, the plan form is unyawed, Jf?= 
4£/7r and 

I) _ 1 _ 1 

D k ~ (1 +0 2 £ 2 ) ^“[1+ (wpdR/4 :) 2 ] * » 

as given in equation (27). This relation furnishes the values 
of DID k in figure 10 at ^ = 90°. If 5£=1 and £^0 one has 
J?= °o, ^=arc tan £ and 

D_ = _ 1 

D K { 1— 0 2 tan V)^ 

This is the general drag relation for the yawed wing when 
aspect ratio becomes infinite. It is to be noted that drag 
remains finite except when the angle of yaw is equal to the 
free-stream Mach angle. In figure 10 the drag curve for 
infinite aspect ratio must therefore have a singularity at 
^=45°. If 5£-> 1 and £-^0 so that aspect ratio remains 
finite, it follows that >0 and D/D K = 1. 

General results in connection with minimum drag with 
fixed base area. — The connection between the Lagrange 
multipliers of the variational problems (see eqs. (12) through 
(15)) and certain quantities in the combined flow field has 
been mentioned. Another interesting interpretation for the 
parameter X can be found. In the case when /= const, on 
Ci, X can be evaluated explicitly in terms of the geometry in 
the plane of T from equations (12b) and (22) 


tan 2\p= (39) 

Figure 10 is a plot of DID K against angle of yaw for 
M= a 2 and J?=4/7r, 4, 8. The smallest of these values 
of aspect ratio corresponds to a circular plan form and ob- 
viously must be independent of the drag ratio is D/D K = 
a / 2/2 and this is in agreement with equation (27) for the 
special case of the circular wing. Several limiting forms of 
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Figure 10. — Optimum drag for yawed elliptic wings. 


X 

(Zoo 



(40) 


Equation (40) states that — X/g*, is equal to the pressure co- 
efficient predicted on the linearized theory for the unidimen- 
sional duct flow bounded internally and externally by the 
characteristic traces, that is, by the curves C x and C 2 , as in 
figure 2. By the other interpretation, \/q m is also the pressure 
coefficient in the combined flow field. This result can also 
be generalized to the case where x=f=my-\-b 0 bv sweeping 
back the entrance to the duct by the same amount as the 
plane of IV 

Another point that should be made in regard to the 
optimum configurations for which /= 0 on T x follows from 
equation (22). That is, the Karman ogive has the greatest 
value of minimum drag for given length, base area, and Mach 
number. This follows directly from equation (22) by noting 
that the area S is a minimum when the trace C x is the circle 
connected with the Karman ogive of the prescribed length 
and base area. The curves in figure 8 show this result clearly. 

The final result of a general nature to be noted here is 
found by comparing equations (12a) and (14a), the mini- 
mizing conditions for given base area and for given lift, re- 
spectively. If differences of over-all sign (because of the 
quadratic dependence of drag upon x) are disregarded, the 
two problems will lead to the same x> hence the same drag, if 

V 2 /— 0 in S 

2g 00 <r=lim (X/ 2 ) on C 2 

z -»0 
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It is found that a supersonic-edged wing with a straight 
supersonic trailing edge, parallel to x=ky, leads to a function 
/ which satisfies equations (41). If the trailing edge is 

x—ky=a (k<CP) 

then the surface 2 2 (see fig. 1) is an inclined plane: 
x =f(y,z) = ky — (iff 2 — k 2 ) H z+a (z> 0 ) 

so that 

/.“-Os*-**)* 

v 2 /=o 

From the second of equations (41), 

2?»<r=-(/J 2 -F) w X 

The drags of the two cases are equal, giving, by equations 
(14b) and (12b), 

(i 3 2 -k 2 )» L 

2 (Zoo 


A=- 


Thus one has the result that if the minimum drag due to lift 
of a supersonic-edged wing with a straight supersonic trailing 
edge parallel to x=ky is known, then a nonlifting wing of the 
same plan form and with base area given by 

A _ ( P 2 —k 2 ) K> L 
2 q~ 

will have the same optimum drag. 

WINGS ANI) BODIES WITH GIVEN VOLUME 

The variational problem in terms of x for the case of 
prescribed volume with zero base area (eqs. (13)) contains 
an arbitrary constant k whose magnitude is to be deter- 
mined by application of the zero base-area condition to the 
solution of the boundary-value problem. Thus set 

X+ 2 P ~V [/*+*/+? (y 2 +z 2 )-\=n(y,z) 

The problem in terms of 0 is then, by equations (13), 

V 2 S2= 0 in S') 


bn = Atft* be 

dn 4 p ,r , dll 


on C 2 


(42) 


a= 2 P !crA p+lcf+ ^J onC \ 

The closure condition can be expressed as 

JJ LfyXy+M dy dz = 0 

s' 

which becomes, after an application of Green’s theorem, 

9 kS=— \ ,f~ (<l- Ads (43) 
2 pJJ* Jc, dn\ 4 pJJ" / 

Clearly now, if /= 0 on Ci, also &=0. Again, as in the given 
base-area case, proper positioning of the origin can simplify 
the analvsis. 


Quasi-cylindrical body of given volume. — The body and 
notation are shown in figure 11. The external wave drag is 
to be minimized under the conditions that volume is fixed 
and base area is zero. 





H > "I 

Figure 11. — Optimum body of revolution with given volume. 

In this case, S 2 is the surface of the conical frustum, whose 
equation is 

x =f( r )=l3(h+R-r);h=a W 

Since /=0 on (\ in this case, the remarks on k , just above, 
apply, and the boundary-value problem is 


K* + 2 7jr.rV-* 

0+i 


d 

dn 


x=0 


in $ 

on (\ 
on Ci „ 


(44a) 


Since (\ and C 2 are circles, and the problem is independent 
of the peripheral variable 0, the solution of equations (44a) 
can be written in the form 


nP / 2i/i / 

X+ o r T T z ; V (ar 2 +b+cln ) 

^PoqC oo P co oo \ 1 / 


(44b) 


where a, 6, c are to be determined from the boundary condi- 
tions. One easilv finds that 


a=-\; 6=I(A>+A) 2 ; c=-\lP 


The drag is then determined as 

D =yJJ W+xS)dy dz=T P „J rx r 

S 

Trgy 2 f h(h+ 

8^00 \ 4 


R-j-h 

rx r 2 dr 

R 


= *§vf h(h± . 2R ) [(/,»_)_ A) 2 — 3/? 2 I + 7?4 
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From equation (13b), the parameter ju can be replaced b} 7 
2D /V, giving finally 


D_2V* f h(h+2R) 

\ 4 


[{R+h) 2 -\iR 2 ]+Knn 



U* l 
/■ ( o 


(45) 


is also an ellipse. Figure 12 shows the wing, the boundary 
curves (\ and C\, and the region S. Let C\ be given by 

£+!*=l 

lP'C* 

I- 2 c A 


where a=pR/l, and 

C(ff) ! < 1 +4<r)[(l +2o-) 2 — 1 2cr 2 ] -f 64<r 4 Z n (46) 

The result of equation (45), as in the case of given base 
area, covers the entire spectrum of fineness ratios and yields, 
in its limiting forms, the results of two-dimensional airfoil 
theory (biconvex section) and slender-body theory. The 
latter case, which is the Sears-Haack slender body with drag 
D S - h (refs. 14 and 1 5) , corresponds to (7=0. Equation (45) 
then becomes 

P s - H = 128F 2 

(Zoo Tfl 4 

The above problem was considered previously by Heaslet 
and Fuller (ref. 16) without recourse to the present techniques 
but, rather, by minimizing after expressing drag in terms of 
the source distribution that could be assumed to generate 
the external shape of the body. In this approach, it becomes 
necessary to find first the source-distribution function, under 
minimizing conditions, and to calculate drag and volume 
subsequently. The details of the calculation are thus less 
direct since the desired quantities are expressed as integrals 
involving the hyperbolic influence function of the supersonic 
source. In reference 16, the function C(a) of equation (46) 
appeared in the form (in a modified notation) 

C(a) = ]^ [( t 7 + 2 < 7)(1 — rj+2a)] ] -[r}(\— rj)E— 

<r( 1 — 4(7) (K — E)](Itj 



2 B H 

Figure 12. — Elliptic wing and characteristic trace. 

The curve C 2 is now merely a segment of the y axis. The 
solution of the problem where minimum drag is sought, for 
given volume and zero base area, is eased by the fact that 
by manipulation of equation (6b), using the extremum 
conditions of equations (13) (with k = 0), the expression 
for the volume can be put in the form 


where K and E are elliptic integrals of first and second kind, 
respectively, of modulus 

/ T 17(1—17) T* 

L(r? + 2(r)(l — >7+2(7) J 

The immediate advantage of equation (46) is, of course, the 
natural one provided by any analytic representation with 
its precise determination of magnitude and rate of change. 
From a disparate point of view, the equivalence of the two 
results gives not only a new fundamental identity in the 
theory of elliptic functions, but also indicates a method 
whereby further identities may be generated. From the 
standpoint of direct application, however, the results of 
reference 16 remain unmodified. The calculations that were 
used to plot the variation of C{o) were found to check to at 
least four significant figures with the present formula, and 
thus provided a satisfying confirmation of the numerical 
techniques used in the original evaluation. 

Elliptic wings of given volume. -As noted previously, a 
wing of elliptic plan form leads to a boundary curve C\ that 



s 


Thus, explicit determination of x from equations (Fla) is 
avoided, and the complicated expression for f need not be 
exhibited; only the fact that it vanishes on the outer bound- 
ary curve (\ is required. It is found that 



fi fPBPC 2 / y 2 z 2 \ 

2 V 


(48) 


satisfies equations (13a) (with lc= 0). Substitution of equa- 
tion (48) in equation (47) yields 



4 


B 2 C 2 S 
B 2 +C 2 


where S ( =ttBC ) is the area enclosed by C\. Now, from 
equation (13b), the drag can be evaluated as 


Z? _4F 2 /J_ 1 \ 

goo pS\B* + C y 


(49a) 
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Finally, in terms of llic original wing parameters (sec 
fig. 12), 

/j*=o*+j8-V 

c 2 =/r 2 c 2 

C? 

R2 4-9 -- 
D 4V 2 P 1 a 2 

(• irac)c 2 (j3 2 +c 2 /a 2 )K ( 49b ) 

which agrees with the result of Jones in reference G. (A 
typographical error in the reference has been corrected in 
equation (40b).) 

It might be supposed that, once again, a simple extension 
of the above results would lead to the minimum-drag result 
for a yawed elliptic wing of given volume. This, unfortu- 
nately, is not the case. The drag found by the present 
method is much lower than the known value, being correct 
only in the limiting case of zero yaw, as above. The reason 
is, presumably, that strict closure along the trailing edge 
lias not been enforced, only the condition that base area 
vanish. An unreal wing has therefore been evolved in the 
yawed case, with patches of negative base area where upper 
and lower surfaces have crossed. It does not seem feasible 
to enforce pointwise closure of the trailing edge in conjunc- 
tion with the present method for drag calculation. Thus, 
caution is necessary when using this method for given volume 
problems. Of course, this difficulty does not arise when 
treating lifting surfaces, for in that case, a negative ordinate 
for the upper surface is of no concern. 

CASES OF GIVEN LIFT ANI) MOMENT 




Figure 13. — Hyperbolic-edged wing and characteristic trace. 



V 2 (x+U a Tzf)=0 

in S ' 

-1 

bn | 

x + U „ {a — y m t) 2 + U „ T2/J = 0 

on C 2 ► 


X=0 

on <7, _ 


(51) 


where x m is the coordinate of the center of pressure. A 
simple, exact solution of equations (51) follows directly if 


A family of wings with supersonic edges. -Consider the 
family of wings whose plan forms are all portions of the 
hyperbola asymptotic to Mach lines through the point 
(— rf,0,0). The equation of the leading edge is 

j} 2 y 2 =2dx-\-x 2 

and the trailing edge is 

x=l 


in which case 




1 U„pL z_ 
162c [l(2d+l)¥ 


W+$ d+l) (20|*| 


/)] 


(52) 


In this event, the drag is given by 

d+jl 

< l <* 4 p \qj [1(2 d+l)]» 


(53) 


where the quantities d and l are shown in figure 13. The 
root chord of the resulting wing is l. If d— >0, the wing be- 
comes a triangle with sonic edges, and if d/l» 1, the wing 
has very large span compared with its chord. 

The surface (fig. 1) is composed of a pair of inclined 
planes 

Xz =f(y,z) =l—p M (50) a 

and the boundary curve C x (fig. 13) is made up of two 
parabolas 

0V=(2rf+/)(/-20|* |) (50b) 

If minimum drag for fixed lift and center of pressure is 
sought for the wings of this family, the variation leads to 
the problem 


where L is the given lift. Since the wing area is 
t | \d+l) jl(2d + l)-(P COsh -1 
the drag parameter is, written in terms of Z=d/l 

/fV 2 = 40 (1+^-)** + Vl+2f— f 2 cosh -1 (1 + f)] (54) 

Figure 14 shows the variation of the drag with f. This 
latter parameter is, in geometrical terms, 0%//, where p 0 is 
the radius of curvature of the leading edge at the apex, and 
l is root chord. As f varies from 0 to <», the plan form 
ranges from a sonic-edged triangle to a wing that has nearly 
a parabolic leading edge. 
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Iii the limit f— >0, when a sonic-edged triangular wing re- 
sults, the value of the drag parameter given by equation 
(54) is 


Cj> 

PC L * 


=0.225 


(55) 


This value is in agreement with the result of reference 8 for 
center of pressure at 60 percent of root chord. Also, the 
approximate result for given lift alone from reference 10 is 
quite close to that of equation (55), being 0.223. This 
would indicate that the center of pressure is near the 60- 
percent-chord position for given lift. In fact, from the 
results of reference 8, it is found that the center of pressure 
for given lift alone lies at 63 percent of the chord. 

At the other extreme, the wing plan form is very 

nearly a parabolic segment. The center of area of such a 
parabolic segment lies at 60 percent of the chord, which 
indicates that the loading is uniform over the wing. This 
is the correct result for minimum drag with given lift in 
two-dimensional flow, and the wing is indeed becoming 
nearly a two-dimensional case at f— ><» . 


DETERMINATION OF SURFACE SHAPE OF OPTIMUM BODIES OF 
REVOLUTION 

Generally speaking, the information gained by solving a 
minimum-drag problem according to the method used in 
this report includes knowledge of the function x> which is 
of course the perturbation potential evaluated in the 
rear enveloping characteristic surface 2 2 of figure 1. It 
would lie useful if this information could be used then to 
determine the singularity distribution, and ultimately the 
surface shape, that gives rise to the flow field having the 
optimal properties in question. In general, however, no 
such method at present exists, 2 although various approaches 
involving some degree of approximation have been indicated 
(ref. 8). In one case, where the singularities lie on a single 
line, a complete analytical solution is possible, and will be 
considered next. 

Ducted body of revolution with prescribed base area. The 
potential due to a line distribution of sources is given by 



tribution; B{x) — 0 for x<—fiR. In this, R is the radius of 
the control cylinder on which the boundary conditions are 
satisfied; the initial point of the body is z=0 when the sources 
start at x=—@R. 

In this coordinate system, the surface S 2 is the conical 
frustum given by 


x =f( r ) = l + l3R-p r =p(2h+R-r); 2 0h=l 

The function x is found by inserting this value for x in the 
above potential relation; 



J /3(2A+«-2r) 
-0/2 


B(xi)dxi 

\0(2h+R)-ztf* b8(2A+ JB— 2r) — 


Now set 


(56a) 


B(z l ) = [p(2h+R)— x x \ l/i Bud 


p(2h+R)-2Pr=t; 


P(2h+R)-t * 

0 


(56 b) 


(r i r~ 0r rU'Qdjt , 

* 27 rj_ w Kx-xJ*- 

where R(x) is the source strength per unit length, and the 
lower limit, x=—/3R, indicates the starting point for the dis- 


2 Since the completion of the work of this report, the authors have seen a paper by E. W. 
Graham, “A Geometrical Problem Related to the Optimum Distribution of Lift on a Planar 
Wing in Supersonic Flow,” Rep. SM-23020, Douglas Aircraft Co., Nov. 1957. In this paper, 
a solution is given to the problem of determining lift distribution on a wing when oblique 
line integrals of the loading taken across the wing are known. Now the lines along which 
such integrals are known all lie in the supersonic “zone of silence” associated with a given 
point, and Graham makes a continuation of the function representing the integrated loading 
as a function of the angle of the oblique line so that the integrated loading can be considered 
known for all lines through a given point. The integral equation expressing the integrated 
loading in terms of local loading is then inverted. Presumably, an analogous procedure will 
hold for the thickness case. For bodies of revolution, however, all oblique lines cut the singu- 
larity distribution in a point, no dependence on the angle of obliquity exists, and the inverse 
problem of finding the singularity distribution can be solved without difficulty. 


and equation (56a) becomes 

_J_ f* B(xi)dxi _ _J_ C r 

(t—Xt)^ 2ir Jo (r—x 2 )^ 


where and x 2 =X\-\-flR. Further, since ^=—2/3 

Cjv * r 

the last equation can be written 


yMti pB U),u 2 

1 tMt Jo {r—X 2 ) ’ A 


(57) 


Equation (57) is recognized as the dual relation to the Abel 
integral equation, namely 



X, (x 2 ) dx 2 
(r-*,)» 


( 58 ) 
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Thus, when Xr is known, the evaluation of B (x), hence of 
B(x), is immediate. 

In order to determine x and hence Xr , one has, from equa- 
tions (12a) 


v '( x+ s db 

f)= 0 

in 7?<r<A > +A 

d / , X 

/)=o 


c>,( X+ 2 P „f7. 

at r= /? 

• 

X=0 

at r=I\-\-h 


Thus, the solution for x is readily seen to be 
X= 9 x ^. (r-R-h) 

IPcoUoo 

and so 

X/9 

Y = — 

r 2 p U 

^ r oo'-' oo 

The solution for t lie source-distribution function can now be 
found by using equation (58); 

B ( r) = — jy- Vt= — fr- V t+PR 

P OO t' 00 P 00 t 00 

Then 

B (x)=-j r [(x+pR) (l+pR-x]» 

which checks with the result of reference 12, since, by equa- 
tions (12b) and (21), it is found that 

y PjJJA 

p 2 S 

and so 

B(x)= U ^[(x+pR) ( 1+pR-x )P (59) 

The surface shape of the body is now readily found, and this 
has been done in reference 16. 

Ducted body of revolution of given volume. —For this 
case, it is convenient to place the origin at the center of the 
body, since x is already known from equation (44b). The 
same solution, as given by equation (58) with appropriate 
changes, can be used, and the resulting source-distribution 
function is 


B [(x+pR) ( l+pR-x))* - 

2 ^ !cos '‘/t3s} 

where for this expression the source distribution starts at 
x=—pR as in the last example. From equations (13b) and 
(45), one finds that 


2 Pm U n 


U V 
2l 4 C(<r) 


where C (<r) is defined in equation (46). The source strength 
is then 


B{X)= 2 

2 ^ ,C03 "rrm} (m> 

Again, determination of the shape is made in reference 16. 

DETERMINATION OF OBLIQUE INTEGRATED LOADINGS 

A procedure has been outlined for relating integrals of the 
perturbation potential <p taken across various lines in the 
surface of a wing and in the rear characteristic surface S 2 
(see, e. g., eq. (19b)). The analysis will now be carried out 
in detail for the family of hyperbolic-edged wings treated 
above. For the special case of the sonic-edged member of the 
family, some of the geometrical relationships are shown in 
figure 3. The auxiliary plane S 3 (given by eq. (17)) is seen 
to cut the wing plane 2 0 in the line CD, and to cut the lower 
part of the S 2 surface in the line EF . Clearly, as the inclina- 
tion of CD, or the chordwise position of # 0 , is changed, the 
line will eventually cease cutting the right leading edge of 
the wing and intersect the wing trailing edge. These regimes 
are distinguished by dividing the range of r 0 into Region 1 
and Region II (as shown in fig. 15), whose bounds are 
functions of the angle 6; 



Figure 15. — Dividing points distinguishing regions I and II. 
Region I : A^ 0 = — d ( 1 — sin 6) < x 0 < X x = l — V l 2 +2dl cos 6 

K61) 

Region IT : X l <x 0 <X 2 =l+y/T 2 ^Ml cos 9 J 

where d and l are defined as in figure 13, and cos 0=f3 tan /i 0 . 
The regions will be considered in order. 

Region I. — The integral relation that results from applica- 
tion of Green’s theorem to the region bounded by 2 3 and 
2 0 (see fig 3) and for which z ^> 0 has been given as equation 
(19b). If Green’s theorem (eq. (18a)) is applied to the lower 
portion, the following relation results: 

-Ijf^dy+fj^dS+ff^dZ-O 

2o 2*3 £2 


483644—50 
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Using again the coordinates v ) s in the separate planes, one 
finds 

J > *»(2 o 2,)<i«J+J%(2 1 Z,)d«=0 (62) 


line in question lying somewhere in the “zone of silence” 
corresponding to the point (a- 0 , 0 , 0 ). Once again, these 
results are not confined to the case of optimum wings. 

Now return to equations (64). It is necessary to determine 
the argument of v in the last integral on the right in equation 
(64b), and this is done by solving for x and z in terms of y 
from the equations of 2 2 and 2 3 ; 


The notation ip (2 2 2 3 ), etc., means that ip is to be evaluated 
along the intersection of the surfaces 2 2 and 2 3 . 

In order to express the element ds in terms of, say, dy 
along the various lines, one can set up the coordinate trans- 
formation to express the N,v,s system as a rigid rotation of 
the x,y, z system. It is thereby found that along CD , ds=dy 
sin d, and along KF, ds=dy. The relations (19b) and (62) 
then become 


J'J <p//jv/?/+sin0j <p(2o2 3 )dy=0 (s> 0) 


(63a) 


-1 j'j<pjlxdy+2{y(2 2 2 3 )dy-smdjy(X 0 2 3 )<ly=0 (z<0) 

(63b) 

where the 2 ^> 0 , z<C0 are necessary because is discont inuous 
at the z=0 plane. 

If the limits of integration for the surface integral over 
2 0 are inserted and the equations differentiated with respect 
to Xq (recall that <p=0 on the surface 20, equations (63) 
become 


One finds that 


2 2 : x=l+/3z 

2 3 : x=x 0 -\~P cos dy — ft sin dz , 


^ ^ x _ R O+Xo+P cos dy — l-\~x 0 -\~P cos dy~ \ 

<^(2 2 2 3 ) — sin e -,y, p (l+sin$) J 

If the differentiation with respect to x 0 is now performed in 
equation (64b), and the result subtracted from equation 
(64a), one gets t lie relation 


j; 


c » (ab+g cos e y>y^y = — p sin gj+sm d “ 


J (0u-\-w)dy 




(3 sin 0(l+sin d) 


(67) 


Making use of the definition of integrated loading given in 
equation (65), one has finally 

L(x O) 0)_ 4 


<U 


U m p sin 0(l+sin d) 


f, x4v 


(68) 


J <Pz(* o+P cos 6y,y)dy = —0 sin 0 | <p z (x n -\-0 cos 6y,y)dy 


(64a) 


( „+0 cos Oy,y)dy=—0 sin 6 | <p x (r 0 +p cos 6y,y)dy+ 

Jc JC 

20 | ' y [^(2 2 S 3 )dy] (64b) 

J E 


where once again the results apply respectively to the cases 
where z approaches zero from above and from below. 

Equation (64a) is seen to be a relation between integrals 
of perturbation components taken along the same oblique 
line across a supersonic-edged wing. If the integrated 
loading 


LM) 1 

(Zco u c 



(65) 


is introduced, equation (64a) becomes 3 


L(x 0 ,d) 

<L*> 


U m p sin 0, 


j w(jr a +0 cos 6y,y)dy 


(66) 


where w is the vertical perturbation component. Thus, the 
integrated loading is proportional to the averaged wing 
slope, both taken along the same line across the wing; the 

3 The authors are indebted to Professor P. A. Lagerstroin and Dr. M. E. Graham for point- 
ing out an error in the definition of the line element ds which introduced an incorrect depend- 
ence on 6 in the versions of the present formulas (66), (68), and (70) given in the original 
NACA Technical Note 4227 (see footnote 1). 


which relates the integrated loading for a point in region I 
to an integral of the function Xz across the rear enveloping 
surface 2 2 . This form of the result is well adapted to the 
case of optimum wings since, in general, surface values of 
the slope are not known, whereas, in the present method of 
optimizing, the function x is known. 

Region II. The plane 2 3 now intersects the trailing edge 
of the wing, and the boundary of the region of integration 
will include a portion of the rear enveloping surface Z 2 for 
both z ^> 0 and 2 < 0 . Figure 16 shows a trace in an £=const. 
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plane of the lines of intersection of the cutting plane 2 3 with 
both upper and lower portions of S 2 . The oblique line across 
the wing intersecting the trailing edge is also shown, and 
the figure serves to define the notation used in this section. 

The formulas resulting from application of equation (18a) 
to both upper and lower portions of the volume bounded by 
the forward Mach surface, the rear enveloping surface, and 
the plane 2 3 are as follows: 

J JJ <pJxdy+ 2 S 3 )cfe+ 

*S 0 

(%&***+ r%(so2 2 )&= 

Jc J D 

JJ <Pzdx dy+2f%(2 a 2 t )ds- 

f ^(2 0 2 3 )&+ f%(2 0 2 2 )cfe= o 

Jc Jd 


=0 


(z> 0) (69a) 


(z< 0) (69b) 


In these equations, the limit indicated by an asterisk is the 
intersection of leading and trailing edges of the wing. Once 
again the equations are differentiated with respect to x 0 and 
the results combined. Thus one finds for a point x 0 in 
region II (see eqs. (61)) 


L(x 0 ,6) ^ 4 

q o> (Jocose 


(j^e’ °) + 


/ f X/ly X z dy \ 

\1— sin0 — i+sin0/ ^ 70) 


U sin 0 \ 1 — si 

For the case of interest here, Xz is a constant, as is seen 
from equation (52). After substituting the values for x and 
Xz in equations (68) and (70), and using the dimensionless 
variables 


x 0 By d 

~T V ~T’ f= 7 


one gets the following results for the integrated loading: 
Region I : 

lUH o,0) 3 


Vf Ve 


(71a) 


Region II: 

ILjj o,») 


3 


L 


4(l+2f) 


n 


2yl l+2f sin 0(l+sin 0) 

(1 +2f)— t?z > 2 , 3 1— sin 0 T l+sin 6 


Vh V D _|_ V D V E 


COS 0 


2Vl+2f 


sin 6 


(71b) 


where 

ij F =(l+sin 0)- , {(l+2f) cos 6+ 

V'2(l + 2f)(l+sin 0)[f o +f(l — sin 6)]} 
i»*=(l+sin 0) -1 {(l+2f) cos 0— 

- v , 2(l+2f)(l+sin 0) [£ 0 +f(l- 


-sin 


«)]} 

r? w =(l— sin 0)~ l {(l+2f) cos 0— 

V2(l+2f)(l— sin 0)[£o+7(l+sin 0)]} 

1?D=(1— &)/C0S 0 


In figure 17 these integrated loadings are drawn for the 
three members of the family corresponding to £=d/l= 0, 
1.5, 5, with 13= 1 and ju 0 = 0°, 15°, 30°, 45° in each case. 
The dimensionless quantity plotted is 

P(£o,Q)=IL(£o,0)/L 



C 1.0 2.0 3.0 4.0 

to 


Figure 17. — Integrated loadings for family of wings with hyperbolic 

edges. 
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It is to be observed that when /z = 0°, 0=tI2 and region II 
ceases to exist. In this case, the integrated loading is 
precisely the chord loading of the wing. One then gets 

IL&,0) 3 t 
L = 2 n ^ () 

independent of the shape parameter f. 

Ames Aeronautical Laboratory 

National Advisory Committee for Aeronautics 
Moffett Field, Calif., Nov. 20, 1057 
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